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1. Introduction 

Recent evidence is piling up in favour of a dark energy component in 
the dynamics of the universe forming up to 70% of the energy den- 
sity: from SNIa magnitudes [Riess et al. 19981 IPerlmutter et oi.~9 98j: and 
more recently |Wang fc Mu kherjee 2004, Nesseri s fc Perivolaropoulos 20041 



Rie ss et al. 2004L |Daly & Djo rgovsk r2004[ lBiesiad a~et al. 2005] : from the 



CMB (Cosmic Microwave Background) measured by WMAP (Wilkinson Mi- 
crowave Anisotropy Probe) coupled with complementary inputs [Bennett et al. 2 003, 
Page et al. 20031 (for reviews, see jPadmanabhan 20031 ISahni 2004] ). The nature of 



this dark energy is still unknown but could lie within two major families: it is either 
coming from a static cosmological constant or developing from a dynamical scalar field 
just emerging from subdominance called quintessence |Ratra & Peebles 1988j . The in- 
terest of quintessence lies in the possibility it holds to solve the four major problems 
raised by the cosmological constant model: the problem of a need for fine-tuning the 
cosmological constant to unnatural values of vacuum energy, the problem of coincidence 
between its present energy density and the order of magnitude of the critical density, the 
problem of the equation of state of the cosmological fluid, which can only take one value 
for a cosmological constant, and the problem of building a model that would naturally 
come up with solutions to these problems while coming from fundamental physics. 

The drive that lead to quintessence models was based on the attempt to address 
those four problems of the observed behaviour of dark energy. The first problem requires 
a theory which yields the present value for the dark energy without the need to 
tune initial conditions at unnaturally small values compared with the natural 
energy scales of the early universe. This addresses the huge discrepancy (typically 
MpJ ' Pde ~ 10 120 ; Mpf. Planck mass; Pde- Dark Energy density) commonly 
found between classical quantum vacuum energy calculations and its cosmological 
measurements. The coincidence problem deals with the current transition status 
of the universe between dominations by matter and dark energy. It requires a 
theory that yields such a tuned peculiar equilibrium without strong constraints on 
its initial conditions, it deals with the value of the dark energy component P/p 

Wang et al. 2"000 



ratio, which is constrained to lie within 



-1 : -0.6 



or even 



within 



-1 ; -0.8 



jEfstathiou 20001 IHannestad ^ M5rtse lT2004j (for dark energy 
equation of state supposed to respect the weak energy condition, thus P/p > — 1). 
Finally the model building problem express the need to have a theoretically motivated 
dark energy potential. 

The fine tuning problem has been the main drive for quintessence models, 
hence it is solved by almost all proposed quintessence potentials, involving high 
energy physics natural energy scales. Various models have been proposed but few 
are deeply grounded in high energy physics. The main examples of motivated 
potentials are the original Ratra-Peebles' potentials [R atra &; Peebles 19 88J 
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that have been linked with the context of global SUSY | Binetruy 1999| 
and the so-called SUGRA potentials derived from supergravity arguments 
[Br ax fc Martin 19991 IBrax fc Martin 2000j . The tracking property is attracting 
much interest because it solves the coincidence problem naturally. It defines classes of 
potentials which Klein-Gordon equation admits a time varying solution that acts as an 
attractor for wide ranges of field initial conditions. In Steinhardt et al. 1999j . general 
conditions to obtain a tracking potential are given. 

In front of the wealth of scalar fields available from physics beyond the 
standard model, the need for discrimination between different quintessential models 
proposed in the literature is a strong drive to confront their predictions with 
observable features. Since the advent of the COBE satellite results and the pursuit 
of refined Cosmic Microwave Background Radiation (CMBR) measurements, the 
attention of cosmologists dealing with quintessence had been focused on CMBR 
anisotropy measurement as primary probes for the various quintessence models (e.g. 

Bra x et al. 2000j ). That has foreshadowed the possibility offered by quintessential 
cosmic dynamics to alter in turn the formation of large scale structures, that 
are more readily available to observations. Some of the earlier attempts in 
this direction that have been started out have restricted themselves to linear or 
perturbation theory of structure formation jBenabed fc Bernardeau 2001] . to pseudo- 
quintessence models (approximation of a constant equation of state parameter different 
than the A term; jLokas fc Hoffman 20011 ILokas et al. 20041 IKuhlen et al. 2005] ) or 
have not pushed the envelope beyond the study of the spherical collapse model 
[Maini nTet al. 2003 bl IMota fc van de Bruck 2004L IJNunes fc Mota 2004j or 
the impact on halo concentrations Dolag et al. 2004] IKuhlen et al. 2 005 . 

The aim of this paper is to compute the non-linear mass function of collapsed 
structures in the presence of a quintessence field. One approach is to use direct 
numerical simulation |Klypin et al. 20031 ILinder fc Jenkins 2003] |Dolag et al. 2004] 
IMaccio et al. 20041 IKuhlen et al. 20051 ISolevi et al. 2005] . but this is only practical 
for a few models with well defined parameters and gives little physical insight into 
the influence of the quintessence potential. Our approach takes into account the fully 
dynamical nature of the field using the methodology developed by Press & Schechter 
[Press fc Schechter 1974] hereafter PS]. This method uses a spherically symmetric 
dynamical model to relate the collapse of massive structures to a density threshold in the 
linearly extrapolated density field. In this way, it is possible to apply Gaussian statistics 
to the initial density field in order to count the numbers of collapsed structures above a 
given mass threshold at a particular epoch. Because the method is nearly analytic, it can 
be applied to rapidly expore the parameter space of possible quintessence models and 
inverted to constrain quintessence models using observational data. Some approaches 
using semi-analytical methods have been published in the course of the present work with 
a restricted range of quintessence potentials jMainini et al. 2003a] IMainini et al. 20 03b, 
Sole vi et al. 2004] INunes fc Mota 2004] . some restricting to very indirect observables 

Ma inini et al. 2003bl INunes fc Mota 2004] . None have yet looked at the same time 
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at observables like the mass function with many different potentials and an approach 
aimed at understanding the dominant physics involved, using the full potential of semi- 
analytical methods. 

The paper is laid out as follows. In section (J2J), homogeneous universe models with 
quintessence are explored and familiarity with the physics of quintessence is developed. 
Some results on the behaviour of the universe's radial scale in various such models 
are presented, showing the possibility offered by large scale structures to probe the 
quintessence potentials. In section (JHJ), the collapse of a top hat spherical model for 
a single primordial inhomogeneity is computed. This provides a simple way of linking 
perturbations in the initial density field with collapsed structures. In section (JH), we 
combined this model with information from the statistics of the initial inhomogeneities. 
Using a PS type scheme we obtain predictions for the mass function in the presence of 
various quintessential models, and compare the evolution of different models. Finally, in 
section El we summarise the results, and explore the constraints that future astronomical 
surveys will be able to set on the form of the quintessence potential. 

2. Homogeneous quintessence evolution 

In this section we will review the features in the homogeneous models of dynamical dark 
energy and emphasize that which points toward an effect on structure formation. This 
roadmap will later be useful for interpreting our results. The first step is to 
establish the behaviour of various quintessence models within FLRW-type solutions. 

2.1. The dynamical system 

In a homogeneous model, the system evolution is entirely defined by the evolution of 
its scale factor and its quintessence field with its time derivative. The Einstein's Field 
Equations governing the Friedman Lemaitre Robertson Walker universe and the Klein- 
Gordon Equation for the quintessence homogeneous, time dependent, scalar field can 
be written as evolutions for the scale factor a and the scalar field Q . 
The first Friedman equation of our model can be expressed as 



the dots refer to cosmic time derivatives and ' will refer (e.g. in Eq 12. 21 below) to 
the total derivative with respect to the field Q. Conventionally, subscript refers to 
the present epoch. Qx = Px/Pc represents the density parameter of species 
X, p c is the critical total density required by a flat universe model. We 
introduce here the density factor Q x = px/Pc 'i note Q Xq = Qx - 

Each term corresponds respectively to the contribution in matter (dark and 
luminous: f2 mo ), radiation (neutrinos and photons: fl ro ), cosmo logical constant (^a ), 
either time dependent (Qq = — ) or present epoch (Qq = Pq / p co , as follows) 



quintessence density factor and the corresponding FLRW curvature f2fc , with the 
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definition 

In this work we have assumed throughout a globally flat geometry of the universe 
(Qk — 0) and a fully dynamical dark energy (^Ao — 0)- 

It should be noted that the very non-linear Klein-Gordon equation requires forward 
integration to reach the tracking solution while the boundary conditions in the density 
parameters are set for a backward time integration. 

The natural units for a primordial quantum field like quintessence are Planck units 
(Planck time: tpi). In those units, expressing the characteristic time, which is the 
Hubble time: f jRahcall et al 19991 IFreeriman 2000] t H \ = H ~ 65km ■ s~ l Mpc~ l ~ 
2.10610~ 18 s _1 = (8.65 x 10 m tpi)~ 1 ), shows a discrepancy between the two timescales 
of the problem of more than 60 orders of magnitudes. Since the Planck time should 
govern the evolution of the field, this could have revealed problematic if we were not 
focusing on cosmic evolution of the quintessence field: neglecting the effects very of rapid 
fluctuations, we can observe the evolution of the field as following a regular evolution 
of its energy density if we adopt the Hubble time as characterizing our system. 

We also assumed spatial homogeneity of the field. [Ma et al. 1999[ 
ILokas fc Hoffman 2001] argue for homogeneous Q, given that they found the small- 
est scale of Q fluctuations to be larger than the clusters scale. Nevertheless it has to be 
mentioned that their results hold their validity from their studies in the mildly non- 
linear regime. jMota fc van de Bruck 2001 IN lines fc Mota 20041 IMaccio et al 2004] 
argue that the highly non-linear regime involved might require some quintessence clus- 
tering or coupling to the dark matter. We will concentrate here on effects in non- 
clustering minimally coupled models as a threshold model to ascertain the impact 
of quintessence on non-linear clustering. 

In Eq. fl2.1|) . where we also assume no curvature nor cosmological constant 
contribution, the effect of quintessence is focused in the density factor, involving 
the energy density of the field Q defined with its potential: in our units it reads 
Qq = ^ + ^pr\ following notations defined above and with Q, the 
quintessence field (Hubble) time derivative, V(Q) its potential energy and 
Tio the Hubble constant in Planck mass units. Its dynamics is governed by the 
Klein-Gordon equation in the case of homogeneity: 

Q= -3-Q-H 2 V\Q). (2.2) 
a 

Though the pressure of the field is not involved in its homogeneous evolution, it is 
crucial to the effect of quintessence on non-linear collapse: for a scalar field recall that 
(with cosmic time) Pq = {^\Q 2 — V(Q)j , and together with the density, they define the 
acceleration state of the universe. One can characterise it using the equation of state 

= Pq/pq (2.3) 
(see figure upper panels). In terms of the energy momentum tensor of the field, 
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Table 1. Our choice of various potentials satisfying the .Ste inhardt et al. 1 999 test. 
They are discussed in section 12.21 They are all either of inverse power, gaussian 
or exponential types. ag and A are slope parameters, Aq characterises the 
potential's energy scale and the values of the field Q are in (dimensionless) 
units of Planck mass, k = 8ttG = &irlm 2 Pl is the gravitational coupling (k = 87r 
in units of mpi). 



Name Potential V 

jRatra fc Peebles 1988j 
SUGRA 

j Brax Martif72000| 

Ferreira fc Joyce 1998| Ag e ~ XQ 



Steinhardt et al 19 99] A.%e~Q 



recall that for scalar fields with these definitions of density and pressure, we have 
Pq = T 00 ; T- = PqSj. 

Now, we will discuss a variety of potentials proposed for quintessence models and 
their previously known homogeneous properties, emphasizing their impact on 
matter domination. 

2.2. Explorations with several tracking potentials 

In this section we restrict our choice upon a set of potentials and recall their previously 
studied equation of state fEq !2.3|) and density parameter (fig) homogeneous evolution, 
emphasizing behaviours that can affect the formation of large scale structures. 

We have narrowed our study on such potentials from the literature that 
we found to agree (at least marginally) with the [Steinhardt et al. 1999j general 
conditions to obtain a tracking potential. We therefore selected several forms of 
potential (the list is not exhaustive) [Ratra fc Peebles 19881 IBrax fc Martin 1 999, 
IBrax Martin 20001 |Ferreira fc Joyce 19981 ISteinhardt et al 1999| for the rest of this 



study. The Ratra-Peebles potential [ Ratra &; Peebles 1988|, hereafter RP or 
Ratra-Peebles] , first discussed potential in the literature has retained its interest 
with its more recent discussion within the context of global SUSY Bi netruy 19 99 . 
The [Rr ax fc Martin 2000| hereafter SUGRA] potential has been motivated in 
the framework of supergravity low energy approximation. The simple exponential 
|Ferreira &: Joyce 1998[ hereafter FJ or Ferreira &c Joyce] potential displays 
a generic form for moduli fields from extradimensional theories flat directions while the 
[Ste inhard t et al. 19991 hereafter Steinhardt et al] potential has been proposed 
as an infinite sum of Ratra- Peebles-type potentials. This set of potential has been chosen 
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as a well motivated starting point. It also spans the main types of potentials, power 
laws, gaussian and exponentials and display very different behaviours. 

Using those potentials, spelled out in table ^ we can explore the impact of 
different models on the timescale of the quintessence dominated epoch and the apparent 
strength of matter domination between radiation and quintessence eras. 

We also can compare the observable effects of various potentials at the homogeneous 
level through the equation of state evolution. The equation of state evolution should 
be attained by SNIa measurements which constrains directly the integrated luminosity 
distance evolution, although it relies on SNIa to be good standard candles without 
systematic errors. In this paper, we are even more interested in the density parameters: 
since, on the rough, dark energy with its negative pressure has a freezing effect on 
clustering of dark matter, the differences between potentials in the equivalence epoch 
for matter-quintessence and in the strength of the matter dominance phase are signs 
respectively of differing inhibiting times for structure formation and of overall matter 
clustering activity. 

We thus emphasize those features in the otherwise known homogeneous evolution 
with Ratra-Peebles potentials for the slope «q = 6 and 11, with the SUGRA 
potentials for the slope «q = 6 and 11, Ferreira & Joyce potential for A = 10 $and the 
|Steinhardt et a l. 1999 potential (the choices for parameter values follow the authors). 
Comparison of the density parameters evolutions allows to conclude on the fact that 
we expect a stronger inhibition of structure formation with the Ratra-Peebles potential 
than with the [Steinhardt et al. 19 99 . than with the SUGRA, than with the Ferreira 
& Joyce potentials. Within models (i.e. for Ratra-Peebles and SUGRA), variations 
of their respective evolution points towards the degree of freedom inherent to 
each potential and thus towards the falsifiability of structure formation tests on each 
model. Another remark drawn from the density parameters concerns the scales of the 
structures affected: in the hierarchical CDM scenario, clusters and superclusters scales 
being formed last, we expect them to be most affected by the changes in inhibition 
epochs from the various potentials because those epochs occur during the most recent 
periods and inhibition is expected to act most on them. Eventually it should be 
stressed that the observational constraints being applied nowadays, differences between 
models are expected to increase as we look back in time. It should be noted that we do 
not expect subdominant quintessence to alter the matter radiation equivalence, and in 
relation, the recombination, thus the power spectrum to remain essentially unchanged 
for the purpose of structure formation. 

The lower panels of figure ^ show some changes for the Ratra-Peebles model 
that are much less pronounced in the case of the SUGRA model, and also illustrate, 
together with the FJ/Steinhardt et al. panel, the variations in equivalence epochs 
and strength of matter domination. The upper panels all illustrate the fact that the 

X This choice is historical and does not allow the marginal tracking behaviour of FJ to 
be reached. It allows nevertheless to explore a very different behaviour of the equation 
of state that yields crucial insights (see section 14. 3|) . 
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Figure 1. Calculations using the Ratra-Peebles R atra fc Peebles 1988) and SUGRA 
|Brax fc Martin 2000] models for the shape parameter values aq — 6 and 11, 
the Ferreira & Joyce |Ferreira fc Joyce 1998| and [Steinhardt et al. 1999) models. 
The equations of state loq of all models are definitely not constant over structure 
formation era. The |Ferreira fc Joyce 1998| model, although constant for a while, also 
varies strongly at recent epochs. The density parameters show (i) matter-radiation 
equivalence is unaffected by subdominant quintessence, (ii) changes between models 
are backwards since they are jointed by present density conditions, (iii) models differ 
in matter quintessence equivalence (which late epoch poses the so called coincidence 
problem) and (iv) models differ in the strength of matter dominance. For both 
Ratra-Peebles models, these discrepancy illustrate the variability of this potential and 
relativize its testability. For both SUGRA models, they illustrate the narrow variability 
of this potential and emphasize its testability, although distinguishing between the 
shape parameters may prove difficult. The equations of state shows that the models 
are well accommodating the measurements |Wang et al. 200 , Efsta thiou 2000| , except 
for the Ratra-Peebles models which is their know pitfall. 



equation of state cannot be considered constant during the epoch of structure formation 
(contrary to |Lokas fc Hoffman 2001] ). 

The solving of Eqs. (j2.1l 12. 2|) was effected with a second order Runge-Kutta 
integration method. As previously mentioned, giving initial conditions for each potential 
shown in table ^ is not constrained with input observations. The search for tracking 
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solutions excludes reverse integration methods: the selection of the tracker solution 
which leads to observed quintessence density was obtained with a simple iteration on 
the characteristic energy scale Aq with respect to the target final quintessential density. 

The cosmic time numerical increment is chosen so as to keep a constant logarithmic 
scale factor increment: alt = a|, where a is fixed (we usually take a = 10~ 2 , a value of 
lCT 3 changing the outcome by less than 1%). 

Given the roadmap that we now have on variations between models in the 
influence of quintessence on large scale structure formation. We are going to effectively 
test them in the mass functions of cosmological haloes within quintessential context as 
computed by the PS scheme [Press & Schechter 1974j . First, in order to get the mass 
function, we need to study the spherical collapse model. This is done in the next section. 

3. Spherical collapse in the presence of quintessence 

The influence of quintessence models on structure formation comes from the repulsive 
gravitational effect of its negative pressure on the bulk of spacetime during its 
phase of less than -1/3 equation of state. Collapsing structures are then slowed 
down in their build up by this relative aggravated expansion. An elementary 
model of this phenomenon can be used to get simple quantitative results that can 
then be inputted in a PS-type scheme that will be discussed in section HJ the 
spherical collapse model (pioneered in [Larson 1969[ IPenston 1969[ IGunn & Gott 1972| 
IFillmore fc Goldreich 1984} |Bertschm gcr 1985) and summarised in [Peebles T9 80j) in 
the presence of quintessence holds the key to this exploration from the beginning of 
the collapse phase after recombination down to shell crossing. After shell crossing, 
mass conservation does not allow to follow the system with just its outer shell but a 
prescription can be used for the virialization of the model. We will now describe the 
dynamics of the cosmological spherical collapse. 

3.1. The dynamical system 

The cosmological spherical collapse is embedded into an FLRW-type universe for which 
all the component of energy density are supplemented with a spherical overdense region. 
Birkhoff 's theorem holds the key to the spherical non linear collapse model: any spherical 
region embedded in a spherically symmetric universe behaves shell by shell as a patch of 
FLRW universe with each characteristics modified to match the corresponding average 
inner ones. In this case, any given shell behaves according to the average of its inner 
density. With a flat background this leads to positive curvature inside the overdensity. 
However, because this inner curvature may change with time, one has to be cautious in 
using the Friedmann's equations [Mota fc van de Bruck 2004| . 

For simplicity reasons and because it is a building block of the original PS scheme, 
we will use the spherical top hat model (sphere of constant overdensity). In this case, the 
averages are equal to the local values and every shell reach the center of the overdensity 
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at the same time, when shell crossing occurs. Problematic for the virialization, this 
feature allows one to follow only the evolution of the outermost shell representing the 
whole system. Virialization is then assumed to take place soon after shell crossing and 
a halo of mass given by the extent of the initial overdensity is formed. 

Thus the radius r d = r of the overdensity can be written as a function of its 
initial value r iod = r, as r od {t) = ri od ak >0 (t) = rjs(t), where we note the scale factor 
of the positive curvature patch a k > as a fiducial or rescaled radius of the overdensity 
region s = r/ri. Thus the radius of the overdensity follows Friedman's equation for the 
modified overdensity patch of universe. For this shifted FLRW model, the Einstein's 
Equations yield the evolution (recall k = 8irG, the gravitational coupling): 

^ (3.1) 

s = -4iCC (3P+p) )] ( 3 - 2 ) 

where the sums on energy density and pressure concern each cosmic species inside the 
patch radius. Again note that the curvature k inside the overdensity is not constant in 
general, so that we prefer to use Eq. (j3.2j) . Since we assumed homogeneous behaviour for 
the scalar field (i.e. no clustering), that means that there is no conservation inside the 
spherical patch for the quintessence field [Mota & van de Bruck 2004] and its pressure 
is that for the background universe. 

The acceleration follows the shifted second Friedman equation (Eq. 13.2)1 . which 
contains as well the pressure terms. The density term for the matter follows the same 
pattern as for Eq. 1)3.1 J) : 

kt 4nGr GM 
-P = ^P 



6 3 r 2 ' 

with the total mass M conserved and contained initially inside the spherical patch. 
The pressure terms depend on the respective state equations. Hence the acceleration 
equation reads, with our time and radius units, 



A - ( Vt ro a 4 + ^(3P Q + p Q 



3H, 



a _lW(l + A0 (3 . 3) 



2s 



Initial conditions of the homogeneous evolution are taken after inflation (a = 10~ 30 ), 
the field is taken in a reasonable range allowing for the tracking solution to establish. For 
the collapse evolution, initial time is chosen in the relevant overdensity linear regime (we 
usually take the arbitrary cut at = 10~ 5 ), the overdensity is set in a Hubble flow, that 
is following the general expansion of the universe at the initial onset of the overdensity, 
(the definition of s sets its initial condition to be = 1 at initial time) so 

it = h ~. (3.4) 

In the model of the spherical collapse, the non-linear density at the 
border of the overdense region can be monitored using the Lagrangian mass 
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Overdensities (linear/ non-linear) evolution 

with SUGRA QCDM model 




a 



Figure 2. Linear and non- linear overdensity evolutions for a range of initial 
overdensities illustrating the computation of the 5 Co (a c ) function for the Ratra-Peebles 
model and several values of A, . The non- linear collapse ends with diverging overdensity 
at the corresponding collapse scale a c while from the linear evolution, initial condition 
Ai is made to correspond to its present day linearly extrapolated value S Co . 



conservation. For the calculation of the overdensity, one just follows the 
canonical definition 5 = pm ~ pb then 



Pb 

3 



5=(l + A i )l—\ -1. (3.5) 
\aisj 

The corresponding initial overdensity value thus proceeds from Eq. ()3.5|) at a = a, 
and s = 1 and using Eq. fl3.4R : 

Si = A h (3.6) 
Si = 3(5, + 1) (- - -) =0. (3.7) 

Here we have used a fiducial initial time from which the collapse was evolved. This 
choice is of course arbitrary and in the PS scheme, this arbitrariness is resolved by the 
extrapolating to present time (redshift z = 0) the initial overdensity 5i = Aj using the 
linear evolution theory 

Q 

5 L + 2~8' L = -Hfo mo a- 3 5 L . (3.8) 

This theory is only valid for overdensities in the linear regime (5 1) but the non-linear 
collapse leads to infinite density at a finite time given by a = a c while the linear theory 
allows the density to remain finite up to z = 0. 

The immediate interest of solving this model lies in the possible comparison between 
different collapse times. The relevant quantity extracted was the function 5 C0 = f(a c ), 
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the linearly extrapolated overdensity as a function of collapse epoch expressed in terms of 
the scale factor. Given a potential, we construct its characteristic extrapolated density 
contrast function of collapse scale: each initial overdensity yields both the collapse 
scale and extrapolated linear scale using the non linear and linear collapse respectively 
(see figure EJ. The construction of the function involved scanning a whole range of 
initial overdensities as illustrated in figure El The critical density contrast function of 
collapse scale factor a c is constructed by the linear extrapolation to the present from the 
arbitrarily chosen initial epoch's starting overdensity. Therefore 5 Co yields an object that 
collapses at an epoch given by the value of scale factor a c . Thus 5 CQ is found by evolving 
our initial conditions with Eq. (|3.8j) 's linear theory, and the corresponding collapse scale 
is found by the non-linear collapse of Eq. (j3.3|) - given by the vertical asymptote at a c , 
in the homogeneous background provided by Eqs. (j2.1| 12. 2|) . 

For implementing this scheme, we first determine the correct QCDM potential 
energy scale, as described in section POl then we bring the model to the field's tracking 
regime and produce initial parameters for the spherical collapse. Eventually the coupled 
evolution of the background FLRW, the quintessence field, the non-linear spherical 
collapse and the linear overdensity evolution obtains the values of S CQ (a c ). 

A fourth order Runge-Kutta integration method was implemented over the whole 
Eqs. fj2.ll 12.21 13.31 13. 8j) system. We still used the logarithmic increment dt = a|, 
where we usually take a = 10~ 2 , as described in section I2~2l but we had to limit the 
lower increment value of dt to avoid inflation of numerical expenses when integration 
approaches the overdensity turnaround point. 

We are now ready to use the Top Hat Spherical collapse model to decide when 
objects are considered to have collapsed, that is when their non-linear overdensity 
diverges. The following section discuss those first results. 

3.2. Critical densities 

Once we have computed a series of S Co (a c ), for a range of initial conditions, set in the 
linear regime, at our arbitrary starting epoch, we can display for each model their 
extrapolated critical overdensity as a function of collapse scale factor and compare them 
among models, which is done in figures E3 (comparing pseudo-quintessence models) and 
|U(for dynamical quintessence). 

Since we are dealing with density contrasts over the background cosmological matter 
density, the main evolution effects are expected to come from the lower homogeneous 
matter density yielded by the quintessence models during their matter dominated era. 
Also the observation setting the models corresponding to z=0 (a = = 1; present 
epoch), we expect the differentiation between models to increase backwards in time. 

A comparison between our various potentials, the ACDM model§ and three pseudo- 

§ It should be noted that for a standard CDM (fi m = 1) model, the spherical collapse yields a linear 
evolution S = & a /a, Peeb les 1980] so, extrapolating the collapse value nowadays, 5 Co = 5 c l/a c which 
gives the straight line asymptote to each models in these log-log coordinates and the well known 
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Figure 3. Comparison between quintessence and pseudo-quintessence models through 
their influence on non-linear collapse scale, function for the linearly extrapolated 
overdensity. The notations are following table ^ with the number referring to the 
index value used in the power of the potential (e.g. RP11: Ratra-Peebles with 
«q = 11). Pseudoquintessence models are referred to by the value of their ad hoc 
equation of state. The spread of curves shows possible distinctions between models. 
The crossing of pseudo-quintessence models by SUGRAll shows none can mimic its 
evolution. LCDM, or ACDM, the boundary cosmological constant model, is shown for 
reference. 

quintessential models (i.e. with constant equation of state) have been performed. The 
overdensity as a function of collapse shows clear distinctions between models are possible. 
The results obtained are plotted on figures EJ for the pseudo quintessences and two 
common models (Ratra-Peebles and SUGRA), and|H for all our studied models. 

The pseudo-quintessence models are shown for reference and comparison to other 
works even though our homogeneous study has pointed that such an assumption is 
highly unlikely to hold compared to real quintessence models (e.g. figure [I}s upper 
panels). However figure El reveals that even though the region scanned by the pseudo- 
quintessence models is covering most of the values taken by our sample of models (except 
for Ratra-Peebles), there is still a strong distinction if measurements are confronted at 
different epoch since no pseudo-quintessence model is matching quintessence evolution; 
for instance, the SUGRAll (notation defined in the figure caption) model cannot be 
fitted with only one wQCDM model. 

The main point is that there are distinctions in this type of representations between 
the various models when examined at earlier epochs and between successive era. Whilst 
evaluation of variability for the SUGRA model yields very little separation, that of 
the Ratra-Peebles are so large that it allows for an observational selection of the free 
power index. Whereas the FJ model crosses the curves from various values of the free 
power index in the Ratra-Peebles model, the |Steinhardt et a l. 1999 model is not very 

S c = 1.686. 




Figure 4. Comparison between quintessence models through their influence on non- 
linear collapse scale as a function of linearly extrapolated overdensity. The notations 
are following table^with the number referring to the free index value used in the power 
of the potential. The spread of curves shows possible distinctions between models (e.g. 
SUGRA6: SUGRA with a Q = 6). LCDM, or ACDM, the boundary cosmological 
constant model, is shown for reference. 



distinguishable from the SUGRA6. However we will see that the sensitivity of the mass 
function to the overdensity allows for some distinction. 

To summarise, it is the departure, in quintessential spherical collapse, from the 
constant value of the collapse critical overdensity of the standard CDM collapse that 
will affect the mass function. On figures El and 01 this is translated into the departure 
from the simple linear log-log relation of the critical overdensity with the scale factor. 
Thus the linear density contrast indexed by the non-linear scale of collapse sets the 
redshift evolution of the PS mass function observable as we will see in the next section. 

4. The mass function in presence of a quintessence field 

In this section we will lay down our assumptions for constructing the mass dispersion 
of our models which imprints initial conditions in the PS scheme, we will then recall it 
before discussing the integrated mass functions for the quintessential models studied. 

4-1. Mass dispersion 

Initial conditions of large scale structure formation are condensed in the PS scheme into 
their mass dispersion. To do this for gaussian fluctuations the basic information lies 
in the density field power spectrum. 

In general, the Fourier power spectrum P(k) for the density fluctuations of scale | 
can be decomposed into a primordial part and a linear evolution part, written as follows 
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f jCole fe Lucchin 1995j pp 267,282, jPeebles 198(1] p169. see discussion in section HX21) : 

P(k, a) = Bk n T Q (k, a) 2 = Bk n T A (k, a) 2 T Q/A (k, a) 2 (4.1) 

where B is a normalization factor and Tq the quintessence transfer function 
written as Tq = T A . Tq/a, T a being the A dark energy transfer function. The 

transfer function is itself written as: T(k, a) =T(k, a = ai).D(a) / D(a,j) where the linear 
growth D gives the time evolution. 

In the following section, we briefly decompose the construction of the structure 
formation initial spectrum between a primordial part and a linear evolution, we 
discuss our assumed spectrum and show the mass dispersion and its normalisation to 
observations. 



4-1.1. primordial power spectrum In the context of inflation, for commonly used 
potentials, one gets a gaussian primordial spectrum P(k,ai n i t .) = \S(k, a in it.)\ 2 =A\ in fik 
(where is the epoch at which a perturbation of size 1=1 /k is reentering inside 
the horizon): i.e the scale invariant Harrison-Zeldovich spectrum, with Xi n fi. and A 
constants depending the inflationary potential and expansion period. 



4-1.2. Linear power spectrum To describe the growth of the fluctuations and the 
evolution of the initial power spectrum P(k, ai n u.) one usually evolves linearly 
fluctuations up to a given epoch a (typically after recombination when - — 1 =z<10 3 ) 
as P(k,t)=P(k, a init ).T 2 (k, a), where T is the above transfer function. 

Previous studies have shown that P(k) depends only weakly on the equation of 
state WQ(a) = P(a)/p(a) |Lokas fc Hoffman 2001j since the Q field is sub-dominant 
at epochs where z is large (e.g., at recombination, z~ 10 3 ). In the case of a 
universe dominated by cold dark matter (CDM) the initial T(k) function is known 
to be well fitted [Bardeen et al. 1986| Sugiyama 1995| . The ACDM (case wq=-1) 
transfer function Ta fit is very close to the CDM function with a slightly smaller 
cutoff scale Efstat hiou et a l. 1992] and identical power laws. Flat cosmologies transfer 
functions fits from N-body simulations of pseudoQuintessence models (!real bf Imodels 
with constant equation of state wq) have been proposed in the form Tq = T A . 
Tq/a |Ma et al. 1999] . They show only modifications to Ta of order unity in the 
observationally relevant range of wq. 

Following previous work, we therefore apply the same initial density fluctuation 
power spectrum to all of the models. This ensures that the differences that we see result 
from the dynamical evolution of the density fluctuations rather than differences in the 
initial power spectrum. In practice, CMBR experiments will establish the fluctuation 
power spectrum independently of the mass functions that we consider here and a 
successful model will have to successfully reproduce both datasets. We thus consider 
that the Ta function given for a model with a cosmological constant can describe suitably 
the fluctuations and do not allow the Q field to influence the initial power spectrum of 
density fluctuations. We therefore adopt the |Bond &: Efstathiou 198 4J form for the 
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power spectrum, following our use of the [Jenki ns""e£ al 20 01] || mass dispertion 
and mass function code, modified with the present work's quintessential 
non-linear collapse program. 

4-1.3. mass dispersion The mass dispersion cr 2 (M, a) is computed using a top-hat filter 
in real space with a radius of filtering R, corresponding to M = 4np b (a)R 3 /3, p&(a) being 
the density of mass of the Universe at a given epoch marked with the scale factor 
a (thus we can express R(M) = (3M/47rp 6 (a)) 1 ^ 3 . 

One has, for the filtered mass dispersion, the relation: 



with the filter Fourier transform W (which is here the top hat filter), indices in Q denotes 
the model dependence. Thus we get the relation: 



a 2 Q (R,a) = a 2 Q {f{a)M l '\a) = D 2 (a) / dk\W(R,k)\ 2 P BondEfsthathlou (k)/87r 3 (4.3) 



with f(a) = (3/47rp;,(a)) 1//3 and D, the linear growing mode for fluctuations. 

4.1.4. spectrum normalization To get the value of the constant in P(k,a) one usually 
normalizes a to the value cr 8 h-i which is observed today for a typical radius R= 8h _1 Mpc 
(here we take h=0.65) and to reproduce the present amplitude of the density contrast 
one gets that cr&h-i is of order 1 (see for example reference |Eke et al. 1996j that yields 

1 /2 

the product cr 8h -i. fl = 0.6, or a sh -i =1.1 for fl = 0.3). 
4-2. The mass function evaluation 

To compute the mass function in the non-linear regime without using fits to N-body 
simulations we restrict here to the popular PS prescription. 

This semi-analytical tool turns out to be the simplest way to construct a mass 
function from the assumed gaussian statistics of the initial density. It neglects a large 
number of effects that tend to compensate each other yielding often accurate fits to the 
numerical simulations. 

It is yet an empirical recipe yielding the number of (collapsed) objects with a mass 
greater than a given one [Press fc Schechter 1974] . 

Hence, we obtain the fraction of collapsed mass linearly extrapolated at a selected 
present time to defined by S Co (t c ) as in section l)3.1|) . and account for Eq. (|4.3j) for the 
variance a itself linearly extrapolated (and normalised with cr 8 ) at the same time to- 
The mass fraction thus writes: 




(4.2) 



F(m > M, a) = 1 - Erf{5 C0 {a)/2a{M, a)), 



(4.4) 



and its derivative with respect to the mass leads to the density of collapsed objects 
|| note their code was built to encompass ACDM dark energy only 
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n PS (m, a) 



- p b (a) (dF(m > M, a) /dM) /M 

(2/tt) 1/2 (p b (a)5 C0 (a)) (da(M, a)/dM) / (Ma 2 (M } a)) 

.exp(-5 2 J(2o- 2 (M,a))) 

(-p b /o-M) (da(M, a)/dM) F PS (M) 



(4.5) 



(4.6) 
(4.7) 



with the characteristic for the PS scheme condensed into the function 



F PS (M) = (2/tt) 1 / 2 (5 co (a)/a(M,a))exp(-5 2 co (a)/(2a 2 (M,a 



))) 



(4.8) 



that can be replaced in later studies with more complex schemes. 

We now have set in place the machinery for computing mass functions for non- 
linear structures in the presence of non-clustering dynamical quintessence. We will then 
present comparisons of integrated mass functions that can be obtained for the potentials 
we selected in the following section. 

4-3. Integrated mass functions 

One of the main interests of this study is to have shown the sensitivity of 5 Co ( a ) as 
function of the chosen Q field potential V(Q). This can be seen above (section ET2"j) with 
figures El and 0] when dealing with the collapse of an initial overdensity. Placing these 
results within the PS scheme, we have found that some differences among the various 
models can emerge from the mass functions. In figure the upper panel shows as a 
reference the mass functions for all the selected models at present time (z=0). The 
models there are all very close to each others, which reflects our z=0 normalisation 
and confirms the usual hypothesis taken of a normalization of the mass functions to 
present days observations. This is why figure |3Js lower panel only recalls the LCDM 
(i.e. ACDM) mass function reference. 

When we examine the mass function at larger z, as shown on figure Efs upper 
panel for z=0.5 and lower panel for z=l, discrimination can occur between the various 
models. All models' mass functions are lying above the LCDM (i.e. ACDM) curve 
for the same epoch. The distinctions remain nevertheless mild at z=0.5, as shown by 
the relative proximity between the curves; they, however, become more pronounced 
at z=l, with the following observations: the FJ model displays the most structure 
suppression, that is the least evolution of mass function, i.e. the least structure 
formation during recent times. Then, in the same decreasing hierarchy, one 
encounters the Ratra-Peebles models, in decreasing order of their power law, the 
Steinhardt et al. model, the SUGRA models, also in decreasing order of their power 
law, before reaching the ACDM model (the Steinhardt et al. and the SUGRA11 models' 
mass functions are almost on top of each other). Surprisingly, those results are not 
completely agreeing with the naive interpretation given in section I2.2[ that would have 
yielded the same hierarchy of structure inhibition and matter dominance except for 
the FJ model, that would be expected to have a lower mass function than all 
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Mass Functions at z=0.5 and 1 



for various QCDM models (and z=0 for Reference) 
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Figure 5. Computation of the integrated mass function for various quintessence 
potentials for past epochs, z=0.5 (a=2/3) and z=l (a=l/2), compared with the present 
epoch (z=0) . The upper panels lack of dispersion among mass functions for all models 
at z=0 illustrates our normalisation and the accuracy of our method. The spreads for 
past epochs, even more so at the earliest epoch, allows for confrontation with cluster 
mass measurements. 
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other quintessence models (see figure lower panels). Confronted with their lookback 
effects on the number of structures present at a certain epoch, these are both agreeing 
with the mass function hierarchy except for the ordering of the FJ model. We are thus 
compelled to consider the only qualitative difference that springs to the eye between 
the FJ model and the rest of them: the evolution of the equation of state (figure HJs 
upper panels). The FJ model appears as the only one with ujq evolving towards higher 
values at its latest stages, during structure formation. It seems thus that the effect of its 
negative pressure in the acceleration equation of the spherical collapse (Eq. 13 .3|) is strong 
enough to counter the strength of matter domination earlier than expected in 
the natural first approach. This can be seen from the ratio of quintessence- 
and matter-induced accelerations of the overdensity, from Eq. (J3.3|) . so we 
get (3 \luq\ — 1)^ > 2 (a + Va) 3 despite having ^ small at a given stage of collapse 
and cosmic evolution. These results cast very interesting light into the intimate 
mechanisms of quintessence that would escape to non semi-analytic approaches. 

The evaluation of model's variability is expected after section !3.2F s more or less 
pronounced differences among each class of model and their respective parameter 
variations (Ratra- Peebles and SUGRA). In the light of integrated mass functions, each 
model evaluated displays little variability so discrimination among power law values 
may prove difficult with this method. 

5. Conclusions 

In this paper we have extended the evidence that several different recent models of 
quintessence predict significant difference in the evolution of the halo mass function, 
in general agreement with previously more restricted works |Mainini et a,l. 20 03a, 
IMainini et al. 2003b(lNunes fc Mota 2004|ISolevi et al. 2004] . and used our wider range 
of potentials to emphasize the impact of the equation of state in the formation 
of dark haloes. Although the best evidence seems to come from the largest mass 
structures, as expected from the homogeneous exploration (section 12. 2|) . the models 
studied do not behave simply according to matter and quintessence dominance 
epochs. Nevertheless, the sensitivity of mass functions to the value of 5 C0 at the 
limits of accuracy cautions us against strong assertions. 

The models can be distinguished if observational measurements of the cluster mass 
function at z ~ 1 achieve a precision of better than ~10% at 10 14 /i -1 M . Because 
the models are normalised to give the same mass function at z = 0, a higher level of 
precision, approaching 1%, is required to distinguish between the models at z = 0.5. 
Several experiments are planned that will exceed this level of accuracy and will have 
good control of the systematic errors. 

One of the most promising is the South Pole Telescope (SPT) survey - a cosmic 
microwave background experiment aimed at detecting the Sunyaev-Zeldovich effect 
(SZE) due to clusters of galaxies [ Ruhl et al. 2004] . The effect arises because inverse 
Compton scattering by the hot intra-cluster plasma distorts the spectrum of the cosmic 
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microwave background. Importantly, the amplitude of the decrement is independent of 
the cluster redshift, and the integrated signal decreases only with the clusters angular 
diameter. The aim of the South Pole Telescope experiment is to map a region of 4000 
deg. 2 , detecting essentially all clusters with masses greater than 2 x 10 14 h~ l M & and 
with redshift less than 2 Carls trom et a,l. 2 002] . The survey is expected to detect 
around 30,000 clusters in the ACDM cosmology, with 30% expected to have redshifts 
greater than 0.8. The amplitude of the SZ decrement is determined by the integrated 
pressure of the intra-cluster plasma, and can thus be used to determine the cluster mass. 
Random errors in the normalisation of the mass function will be ~ 1% at both z ~ 1 
and z ~ 0.5. Systematic errors in the determination of mass and the scatter between 
cluster mass and the observed SZ decrement are a significant source of concern, but 
these can be controlled using the "self calibration" techniques discussed in |Hu 2 003 
and |Majumdar fc Mohr 2003 . An other technique based on the optical detection of 



clusters may also yield promising results when combined with similar self-calibration 
methods |Gladder fc Yee"2005]. 

These observational programmes will map the development of the mass function 
from z = 1 to z = at the level precision that is required to distinguish between the 
quintessence models considered in this paper. The techniques we have developed allow 
this map to be directly related to the form and parameterisation of the quintessence 
potential. 

Although this work focuses on the actual mass functions of studied models, that 
is on the real volume number densities predicted, some authors have argued that 
detection of such density should be hindered by the global geometric effect of dark 
energy [Sol evi et a,l. 2004"| ISolevi et al. 20 05 . However, their treatment of observed 
cluster samples with a unique mass cut-off forgets about the bias-geometry dependence, 
discussed by |Kaiser 1984j and retained as mass-observable relation in dark energy 
studies [Mohr 2004J, that should certainly call for further scrutiny. Such a study 
is in progress that should test whether a break of degeneracy is induced by 
this bias-geometry dependence and that will be at the core of a companion 
paper. 

More potentials have attracted attention |Barreiro et al. 2000] lAlbrecht fc Skordis 2 000 
and although they combine our studied exponential and power laws, they should 
be confronted using our method. Confirmation could be sought using parallel semi- 
analytical methods like the [Jenkins et al 2001] one in the context of quintessence 
[Linder fc Jenkins 2003j . Eventually, the homogeneity and minimal coupling of 
the field lacks proof in the highly non-linear regime IMota fc van de Bruck 2 004, 
INunes fc M ot a 2004 and calls for some wider explorations in the line of 
jMota fc van de Bruck 20041 IJNunes fc Mota 20041 IMaccio et al. 2004j . This will be the 
subject of a follow up paper. 
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